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Discrete quantum walks are dynamical protocols for controlling a single quantum particle. Despite 
of its simplicity, quantum walks display rich topological phenomena and provide one of the simplest 
systems to study and understand topological phases. In this article, we review the physics of discrete 
quantum walks in one and two dimensions in light of topological phenomena and provide elementary 
explanations of topological phases and their physical consequence, namely the existence of boundary 
states. We demonstrate that quantum walks are versatile systems that simulate many topological 
phases whose classifications are known for static Hamiltonians. Furthermore, topological phenomena 
appearing in quantum walks go beyond what has been known in static systems; there are phenomena 
unique to quantum walks, being an example of periodically driven systems, that do not exist in 
static systems. Thus the quantum walks not only provide a powerful tool as a quantum simulator 
for static topological phases but also give unique opportunity to study topological phenomena in 
driven systems. 



I. INTRODUCTION 



Discrete time quantum walk 



Discrete quantum walk, in its simplest form, is a dy- 
namical protocol for controlling a single spin 1 /2 particle 
in one dimensional lattice (see Fig. [I]). It consists of two 
operations given by a spin rotation and spin-dependent 
translation. The evolution of quantum walk results from 
the repeated applications of these two operations in al- 
ternate fashion. From the first introduction of the con- 
cept of quantum walks by Aharono^, quantum walks at- 
tracted tremendous attentions due to their implications 
in quantum information science^. One of the attractive 
features of quantum walks is its simplicity which allows 
any student of physics who has a basic understanding of 
quantum mechanics to grasp its definition. Yet, the con- 
sequence of quantum walks is profound; on one hand, it 
provides a powerful tool for quantum algorithms^! and 
on the other, it displays rich topological phenomena re- 
vealing the deep relation between physics and the ab- 
stract field of mathematics^. In this article, we review 
the topological phenomena appearing in discrete quan- 
tum walks. While the concept of topological phases is 
often challenging to understand because it tends to be 
intimately intertwined with the physics of solid state ma- 
terials, discrete quantum walks provide a rare opportu- 
nity to explain the idea of topological phases in the most 
elementary form due to its simplicity. 

This review article is organized as follows. In Sec- 
tion |nj we study the physics of a one dimensional quan- 
tum walk. First, in Section [II A[ we define the simplest 
one dimensional quantum walk and give the description 
of quantum walks through so-called effective Hamiltoni- 
ans. This quantum walk possesses a symmetry that is 
not apparent at a first sight but plays a crucial role for 
its topological properties, and we describe this symme- 
try in Section II B The effective Hamiltonian approach 



to quantum walks gives an intuition behind the behavior 
of quantum walks, which we use to derive the analytic ex- 
pression for their asymptotic distribution in Section fll C| 
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FIG. 1: The protocol of a conventional discrete quantum walk. 
A conventional quantum walk is a dynamical protocol for con- 
trolling a spin 1/2 quantum particle in one dimensional lat- 
tice. It consists of two operations; 1. rotation around y axis 
by angle 9, whose operator is given by R y (0) = e~ xB,7y ^ 2 \ 2. 
spin-dependent translation T where spin up particle is move 
to the right by one lattice site and spin down particle is moved 
to the left by one lattice site. One step of the quantum walk 
is given by U — TR y (8) and the evolution of the particle after 
many steps are studied. 



This conventional quantum walk and slight variations of 
it has been realized in a numerous experiments with dif- 
ferent physical settings, and we explain some of their 
realizations with cold atoms, photons and ions in Sec- 
tion [udJ 
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In Section [TTTJ we first briefly review the main ideas of 
topological phases that will be relevant to the study of 
quantum walks. 



Starting from Section IV we investigate the topologi- 
cal nature of the quantum walks. In Section [IV A[ we de- 
scribe the concept of topological invariants in the context 
of quantum walks. In order to fully explore the topologi- 
cal phases in quantum walks, we extend the conventional 
quantum walks to so-called split-step quantum walk in 
Section pVB| which displays multiple topological phases 
in its phase diagram. Using split-step quantum walks, 



we argue in Section IV C that physical manifestations 
of topological nature of quantum walks are the appear- 
ance of bound states across distinct topological phases. 
We demonstrate the existence of such bound states in 
two physically distinct situations; one is inhomogeneous 
quantum walks where rotation angles are varied in space 
(Section IV C ) and the other is the quantum walks with 
reflecting boundary (Section IV D I. The unique nature 



of these bound states lies in the robustness of their exis- 
tence against small perturbations. We provide additional 
understandings of this robustness in Section |IVE| which 
is based on the spectrum and topological invariants as- 
sociated with the bound states. These bound states in 
the quantum walk have the same topological origin as 
the bound states predicted to arise in polyacetylene de- 
scribed by the Su-Schrieffer-Heeger model and Jackiw- 
Rebbi model in high energy physics. 

Remarkably, there are also phenomena in quantum 
walks that have not been predicted before in static sys- 
tems; the existence of two flavors of bound states at 
quasi-energy E — and E — ir. This phenomenon is 
unique to periodically driven systems, and we illustrate 
the physics in Section [IV G These one dimensional topo- 
logical phenomena in quantum walks have been experi- 
mentally verified in [5]. 

From Section [Vj we extend the idea of quantum walks 
to two dimension and study their topological properties. 
In Section |V A[ we define the two dimensional quantum 
walks that display non-zero Chern numbers, which is 
the topological invariant responsible for integer quantum 
Hall effects. We explain in detail how Chern number 
arises in quantum walks in this section and demonstrate 
that non-trivial topology in this class results in unidircc- 
tionally propagating modes at the edge of the systems. 
As is the case for one dimensional quantum walks, two di- 
mensional quantum walks display topological phenomena 
unique to periodically driven systems. In Section VB we 



describe a simple quantum walk which possesses unidirec- 
tionally propagating modes as a result of topological in- 
variants unique to driven systems, namely the non-trivial 
winding number in energy direction. Due to its simplic- 
ity, this quantum walk has the advantage of being easier 
to implement in experiments compared to the quantum 
walk introduced in Section IV Al 

In Sect ion \Vl\ we briefly discuss the realization of other 
topological phases, using quantum walks. In Section [VII| 
we conclude with possible extensions of ideas reviewed in 



this article. 



II. ONE DIMENSIONAL QUANTUM WALK 
GENERAL PROPERTY 



A. Effective Hamiltonian description of quantum 
walk 

The simplest form of discrete quantum walks^ is de- 
fined as a protocol acting on a single particle in one di- 
mensional lattice with two internal degrees of freedom. In 
an analogy with spin 1/2 particle, we refer to the internal 
degrees of freedom as as "spin up" and "spin down." This 
quantum walk, which we call as a conventional quantum 
walk in this article, consists of two operations (see Fig. [I]); 

1. rotation of the spin around y axis by angle 9, corre- 
sponding to the operation R y {9) = e~ ierTy / 2 where 
<j y is a Pauli operator. 

2. spin-dependent translation T of the particle, where 
spin up particle is move to the right by one lattice 
site and spin down particle is moved to the left by 
one lattice site.. Explicitly, T = J^j \ j + 1)01 ® I t 

>(t | + |j -i)0|® 11X11. 

These two operations make up one step U = TR y (9) 
of the quantum walk, and the evolution of the particle 
after many steps of the walk is studied. It is possible 
to define more general quantum walks by replacing the 
first operation by any unitary operation R(0,cp), which 
can be written as the product of the rotation of the spin 
around some axis n by angle 9 and the phase accumu- 
lation ip i.e. R{9 1 ip) — e ~ t9a -' J / 2 e - tl P . However, many 
central concepts of topological phases in quantum walks 
can be illustrated with the simple quantum walk defined 
above, so we focus on the conventional quantum walk 
in this paper. The extensions to more general case is 
straightforward. Discrete quantum walks have been real- 
ized in variety of experi ments with ions, cold atoms and 
photons^^. In Section IIP we describe experimental 



realizations of quantum walks in some of these systems 
in details. 

Many properties of quantum walks, such as the dis- 
tribution of the particle after many steps, have been 
exte nsive ly studied from mathematical physics point of 
view^^. In this paper, we take the intuitive picture in 
which quantum walk is considered as a stroboscopic re- 
alization of static effective Hamiltonian, defined through 
the evolution operator of one step quantum walk protocol 
U = TR y {9)- 



U = 



TR y (9) 

c -iH el{ T 



(1) 

(2) 



Here, T is the time it takes to carry out one step of 
the quantum walk. Because n steps of quantum walk 



3 



correspond to the evolution U n = e~ mnT , the evo- 
lution under the quantum walk coincides with the evo- 
lution under the effective Hamiltonian U = TR y (9) at 
integer multiple times of T. From this perspective, quan- 
tum walk provides an unique quantum simulator for the 
static effective Hamiltonian H c g through periodic drive 
of a quantum particle. In the following, we set T = 1. 
While the effective Hamiltonian H c h represents a static 
Hamiltonian, there is one important difference from the 
truly static Hamiltonian; the energy is only defined up 
to 2n/T, as is clear from the definition Eq. This 
periodic structure of the eigenvalues of effective Hamil- 
tonian, called quasi-energy, is the result of discrete time- 
translational symmetry of the system. Just as momen- 
tum becomes quasi-momentum with periodic structure 
when the system possesses only discrete spatial transla- 
tional symmetry, energy also becomes quasi-energy with 
periodic structure in this case. This distinction between 
energy and quasi-energy plays an important role in un- 
derstanding the topological phenomena unique to quan- 
tum walks, or periodically driven systems, as we explain 
later. 

For the quantum walk defined above, the evolution op- 
erator of one step U = TR y {9) possesses spatial trans- 
lational symmetry, and thus the evolution operator be- 
comes diagonal in quasi-momentum space. In particular, 
spin-dependent translation T can be written as 

t = \i + ® I t)(t I + li - 1)01 ® U)U I (3) 



dk e lka * <g> \k)(k\ 



(4) 



In this expression, we see that the spin-dependent trans- 
lation mixes the orbital degrees of freedom represented 
by quasi-momentum k and spin encoded in a z . The pres- 
ence of such spin-orbit coupling is a key to realizing topo- 
logical phases. We note that the continuous quantum 
walk does not have such spin-orbit coupling and thus is 
distinct from discrete quantum walks in the topological 
properties 2 

Now in quasi-momentum space, the effective Hamilto- 
nian for the conventional quantum walk takes the form 



dk [E(k)n{k) -a}(g)\k)(k\, 



(5) 



where a = (a x , <Jy,<J z ) is the vector of Pauli matrices and 
the unit vector n(fc) = (n x ,n y ,n z ) defines the quantiza- 
tion axis for the spinor eigenstates at each momentum k. 
For 9 or 2ir, explicit expressions for E{k) is given by 



cos E(k) = cos(0/2) cosfc 



(6) 



Typical band structure of quasi-energies E(k) is plotted 
in Fig. [2^,) for 9 — tt/2. There are two bands because 
the system has two internal degrees of freedom, and for 
generic values of 9, the two bands are separated by a 
band gap. 
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FIG. 2: a), b) quasi-energy spectrum of the effective Hamil- 
tonian Eq. |5} for conventional quantum walks with rotation 
angle 8 — ty/2 and 8 — 0. The spectrum consists of two bands 
coming from two internal degrees of freedom, and there is a 
finite gap between the two bands for general value of 8 except 
for 8 7^ 0, 2-7T. For 8 — 0, 2n, the spectrum closes the gap as 
is observed in b). Note that the gap is closing in this case at 
quasi-energy E = as well as at E — n. c) The behavior of 
the eigenstates n(fc) in Eq. |7| represented on a Bloch sphere. 
For a given quasi-momentum k, the eigenstate is the superpo- 
sition of spin up and down, and therefore, can be represented 
as a point on a Bloch sphere given by n(k). For a conven- 
tional quantum walk, n(k) traces a circle around the origin 
as k goes from — -n to n. Note that n(fc) is perpendicular to 
a vector A = (cos(0/2), 0, sin(#/2)) for all k in our quantum 
walk. For gapless quantum walk with 8 = 0, 2n, n(k) becomes 
ill-defined at those k corresponding to the gap closing point. 



On the other hand, at 9 — 0, R y {9) = I, and therefore 
the effective Hamiltonian is i? c ff = ka z . Thus the quasi- 
energy bands close the gap at quasi-momentum k = 
at quasi-energy E — 0. Moreover, due to the periodic- 
ity of energy, the spectrum also closes the gap at quasi- 
momentum k = 7r at quasi-energy E = n. The spectrum 
for 9 = is illustrated in Fig. [2]d). Similar situation 
occurs at 9 — 2tt. 

Interesting structure appears in n(fc) of Eq. JH]). For 
9 7^ 0, 2tt, n(k) is given by 



n(fc) 



(sin(6»/2)sinfc, sin(6'/2) cos k, - cos(6»/2) sin k) 
sin E(k) 

(7) 

Note that the eigenstates of the effective Hamiltonian 
H c g(9) are the superposition of spin up and spin down, 
and therefore can be represented as a point on a Bloch 
sphere. The unit vector n(fc) is nothing but the unit 
vector that determines the direction of this point. The 
behavior of n(fc) on a Bloch sphere for 9 = tt/2 as k goes 
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from — 7r to 7r is plotted in Fig. [5J:). We see that n(fe) 
"winds" around the equator of the Bloch sphere. This 
peculiar feature is in fact the origin of topological nature 
of quantum walks, as we see in Section [IV| 

Note that at 8 = 0, 2ir, the states at quasi-momentum 
k = and ir become degenerate and thus any superposi- 
tion of spin up and down becomes the eigenstates of the 
Hamiltonian at these quasi-momentum. Thus, eigenvec- 
tor n(fc) becomes indeterminate at these points. 



B. Hidden symmetry of quantum walks 

The effective Hamiltonian of the quantum walk H e g 
as given in Eq. ([5]) has an interesting symmetry. First, 
we note that n(k) is perpendicular to the vector A = 
(cos(0/2), 0, sin(6>/2)) as is easy to check from Eq. (Q. 
The symmetry of the system is then given by the rotation 
around the axis A by angle n which takes iJ G ff to —H e g, 
or 

r-^effT = -H eS r - 6 i7rA'<r/2 (g) 

Indeed, as is clear from the picture of Fig.[2l such rotation 
takes n(fc) to — n(fe) for each k, and thus takes H c s — 
j dkE(k)n(k) ■ a ® to minus itself. This symmetry 

given by an unitary operator is called sublattice or chiral 
symmetry^. 

One interesting feature of quantum walks which results 
from the sublattice (chiral) symmetry is the symmetric 
spectrum; states with energy E and —E always appear 
in pairs. This is easy to demonstrate. Given a state \ip) 
with eigenvalue E, one can check that the state T\ip) is 
an eigenstate of the Hamiltonian with energy — E: 

H cS \^} = E\4>) 

-► H ee (r\^}) = -E(r\^}) 

There are exceptions to this pairing of states. When E 
satisfies E = —E, the states \tp) and r|i/>) can represent 
the same state. Due to the periodicity of quasi-energy, 
quantum walk has two special energies that satisfies the 
condition E = —E, given by E = and E — it. We will 
later see that this property of E = and E — ir leads 
to the topological protection of E — and E = n states. 
The existence of a single E — state is known in the non- 
driven systems with sublattice (chiral) symmetry-^, but 
E = it energy state is the novel feature of periodically 
driven systems. 



C. Asymptotic distribution; ballistic propagation 

In this section, we provide an intuition behind the 
propagation of a particle under quantum walks by con- 
sidering their asymptotic distributions. Quantum walk 
was originally conceived as quantum analogue of random 
walkPl In fact, it is easy to check that if one carries 



out a measurement of the particle after each step, quan- 
tum walk reduces to a biased random walk. Thus, by 
varying the amount of decoherence in the system, one 
can smoothly go from a quantum walk to classical a ran- 
dom wall^. However, the behavior of a random walk 
is quite different from that of a quantum walk in the 
absence of decoherence. We remind the reader that a 
(non-biased) random walk asymptotically approaches the 
Gaussian distribution with the peak centered around the 
origin where the mean squared travel distance is given by 
(x 2 } — Na 2 where a is the step size of one step. Thus, 
the particle propagates in a diffusive fashion under a clas- 
sical random walk. 

On the other hand, a particl e und er quantum walks 
propagates in a ballistic fasMorPEll. This fact almost 
trivially results from the understanding of quantum walks 
as simulations of static effective Hamiltonian of non- 
interacting particle as given by Eq. ^ and Eq. j5j. If 
one prepares a particle in a state such that its quasi- 
momentum is narrowly concentrated around k, then it is 
intuitively clear that the particle ballistically propagates 
with the group velocity given by v k = where E(k) 

is the quasi-energy given in Eq. ([6]). 

This intuition can be made rigorous by deriving the 
asymptotic distribution of quantum walks. The argu- 
ment above shows that the distribution of the variable 
X = x/N in the asymptotic limit converges to a finite 
form. Here we consider the quantum walk of a parti- 
cle initially localized at the origin x = with spin state 
given by \s). In Appendix [Aj we show that the distribu- 
tion function of X takes the following form 

^{l + {n{k)-a))8{v k -X) 
+ l -{l-{n{k)-a))5{v k +X) (9) 

where (n(fc) • a) — (s\n(k) ■ ajs), and Vf. = 9 gu^ ■ This 
result is quite intuitive; each momentum state k propa- 
gates with velocity ±Vk where + sign is for spin parallel 
to n(fc) and — sign is for spin anti-parallel to n(k). Be- 
cause the initial state is localized at a single site, it is the 
superposition of all quasi-momentum state fc, and thus, 
the asymptotic distribution is given by the sum of the 
contributions for each fc. 

The form of asymptotic distribution written above im- 
mediately leads to various results known in quantum 
walkJSES]. For example, symmetry of the asymptotic dis- 
tribution x O —x exists whenever (n(fc) • a) is an even 
function of k, which is the case when the initial spin state 
is pointing in y direction. Otherwise, the distribution is 
generically not symmetric around the origin. Numerical 
evaluation of asymptotic distribution is always possible, 
and in certain cases, the analytic result can be expressed 
in a compact form. For example, for the asymptotic dis- 
tribution of 9 = 7r/2 quantum walk with initial spin state 
pointing in z direction, the asymptotic distribution is 
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P(X) 




FIG. 3: The asymptotic distribution of a conventional quan- 
tum walk with 9 — n/2, whose analytical solution is obtained 
in Eq. Q and Eq. ( |10[ ). The result is expressed for rescaled 
coordinate X = x/N where N is the total number of quantum 
walk steps. X takes a finite value for N — > oo limit because 
a particle propagates in a ballistic fashion. 



given by 



P{X) = 



1 1 



1 



TTl + X Vl-X 2 



7i- x -T2 (10) 



The distribution is plotted in Fig.[3j The derivation given 
here and Appendix [A] can be easily extended to more 
general quantum walks with different unitary operations 
or even to higher dimensions. 



D. Experimental realizations of quantum walks 

One dimensional quantum walks described in previous 
sections have been realized in experiments. Since quan- 
tum walk is a general concept applicable to many dif- 
ferent physical systems, there are realizations with cold 
atoms, ions, and photon^HUI] Such realizations in differ- 
ent physical settings allow different controls of the sys- 
tems, such as the ability to choose the rotation opera- 
tions, to introduce known amount of dephasing^, or to 
create spatial boundary between regions with different 
rotation operations^. Thus, study of quantum walks in 
experimental settings is versatile and is not usually re- 
stricted by the technology of a specific field. 

The realization of quantum walks takes a widely differ- 
ent forms for different systems. The realization with cold 
atonPis probably the simplest and most straightforward. 
In this experiment, Karsi et al realized the quantum walk, 
using the cesium atoms (Cs). Here, two hyperfine states 
of cesium (Cs) atoms are used as spin degrees of freedom. 
The spin rotation is implemented through the application 



of resonant microwave radiation between these two hy- 
perfine states. The spin-dependent translation, on the 
other hand, is implemented by the adiabatic translation 
of spin dependent optical lattices. This experiment im- 
plemented quantum walks up to 10 steps, and observed 
the distributions of the particle. The experiment shows a 
good agreement, but quantum walks were dephased after 
10 or so steps. 

On the other hand, in the case of photonic realization 
by Broome et aZjlOj . the vertical and horizontal polar- 
ization of a photon is used as spin degrees of freedom. 
The rotation operation is implemented through half wave 
plates, where the polarization of a photon is rotated 
as a photon goes through the plate. The polarization- 
dependent translation is implemented by bircfringent cal- 
cite beam displacer. The optical axis of the calcite prism 
was cut in such way to displace horizontally polarized 
light to the perpendicular direction to the propagation 
direction and transmit the vertical polarized light with- 
out displacement. Now these optical components are put 
in series in, say, z direction, so that the photon goes 
through them one by one as it propagates in z direction. 
Thus, z direction plays the role of time direction, and 
the photon is evolved according to the quantum walk as 
it propagates through these optical components. This 
experiment implemented 6 steps. In the subsequent ex- 
periments, this experiment was extended to create the 
boundary between regions with different rotation angles, 
which was used to investigate the topological nature of 
quantum walks^. Yet another implementations of quan- 
tum walks with photonic architecture is demonstrated in 
[9], using a fiber network loop. 

Lastly, we describe the quantum walk implemented 
with ions. In [7J, Zhringer et al realized quantum walks 
in a phase space with ion 40 c a + They used the inter- 
nal states \Si/2,m = 1/2) and \D 5 / 2 ,m = 3/2) as spin 
degrees of freedom. For spatial degrees of freedom, they 
used the excitation of the ions in the harmonic traps, 
where the superposition of raising and lowering opera- 
tors are identified as the coordinate operator x = a) + a 
and momentum operator p = a + a . In this space, they 
implemented the spin dependent translation by applying 
a bichromatic light that is resonant with both the blue 
and red axial sideband. This shows the interesting fact 
that quantum walks can be encoded in abstract space. 
The experiment demonstrates that the quantum walk in 
this space can maintain the coherence up to even 23 steps. 

While they are not an experimental demonstrations, 
there are two interesting possible realizations of quantum 
walks proposed in natural systems. 

Oka et al proposes in [ISj that the evolution of elec- 
trons on a ring under the application of DC electric 
field can be understood as quantum walks in energy 
level space. When the electric field is understood as 
the time-dependent vector potential, the problem rep- 
resents a time-dependent problem, where the Landau- 
Zener transitions of electrons among different levels can 
be mapped to quantum walks. In this work, they pro- 
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posed the existence of a localized state near the ground 
state of the system, which is a manifestation of the topo- 
logical nature of quantum walks explained in this article 4 
(see Section IV D I. 

Another proposal by Rudner and Levitov-^ concerns 
an extension of quantum walk to include the decaying 
sites at every other lattice sites, which can arise in the 
problem of coupled electron and nuclear spins in quan- 
tum dots in the presence of competing spin-orbit and 
hyperfine interactions. This quantum walk, intriguingly, 
displays topological phenomena as well, where the mean 
walking distance of a particle before it decays at one of 
the sites is quantized to be an integer. 



III. BRIEF INTRODUCTION TO 
TOPOLOGICAL PHASES 

Quantum walks described in Section [n] display rich 
topological phenomena. In this section, we review the 
ideas of topological phases, and provide the background 
for understanding the topological phenomena in quantum 
walks. 

The relation between quantum phases of matter and 
topology was first discovered through the study of integer 
quantum Hall effect, which revealed the quantization of 
Hall conductance for two dimensional electron gas in the 
presence of a strong magnetic fiel d 1 ' ' ls [ The quantiza- 
tion of Hall conductance is very precise, and moreover is 
robust against perturbations such as the impurity of the 
materials. The fundamental origin for such robustness 
is the topological nature of ground state wavefunction. 
Ground state wavefunction of integer quantum Hall sys- 
tems is associated with a topological number. A topo- 
logical number is a global property of the shape of wave- 
function, like "winding," and thus cannot be changed by 
a continuous change. Because Hall conductance is di- 
rectly given by this topological number^, its value does 
not change under the small change of Hamiltonian or its 
ground state wavefunction. 

In addition to the quantized Hall conductance, yet an- 
other direct physical consequence of non-trivial topology 
of wavefunction is the existence of unidirectional edge 
states at the boundary of the sample^. The two phe- 
nomena, the quantized Hall conductance and unidirec- 
tional edge states, are closely related where the current 
for the quantized Hall conductance is carried by the edge 
states. From this point of view, the robustness of the Hall 
conductance against impurities results from the robust- 
ness of unidirectional edge propagation against impurity 
scatterings. As we explain in this article, while quantized 
Hall conductance is a unique phenomenon to the topolog- 
ical class of integer quantum Hall systems, the existence 
of robust edge states is generic feature of any topological 
class. 

Because topology is a general property of ground state 
wavefunction, such idea is extendable to other systems 
in other dimension d 21 l 22 [ In one dimension, Su Schrief- 



fer and Heeger gave a simple model of conducting poly- 
acetylene and found the existence of topological solitons 
at edgeiP^. This so-called Su-Schrieffer-Heeger model of 
polyacetylene is an example of one dimensional topolog- 
ical phase with sublattice (or chiral) symmetry. Inde- 
pendently, physics in the same topological class was also 
studied in the c ontext of high energy physics by Jackiw 
and Rebbi-M25l. j n recen t years, band insulators with 
time-reversal symmetry are predicted to possess topolog- 
ical properties called quantum spin Hall effedpSGH, and 
these so-called topological insulators were soon realized 
in experiments with HgTe 2 ^. One important conceptual 
advance of these topological phases compared to integer 
quantum Hall phase is that these are new topological 
phases appearing in the presence of symmetries such as 
sublattice and time-reversal symmetry, whereas the inte- 
ger quantum Hall effect is the phenomena that appear in 
the absence of these symmetries. The idea has been fur- 
ther extended to three dimensional systems in the pres- 
ence of time-reversal symmetry^HH] 



Motivated by these findings, a several groups inde- 
pendently classified the non-interacting systems in the 
presence of particle- hole, time-reversal, or sublattice (chi- 
ral) symmetry, giving the "periodic table" of topological 
phases^HUl. While the possible existence of topological 
classes are known for the symmetry class within these 
categories, their realizations are not easy in condensed 
mater materials and consequently, some table entry have 
not yet found physical realizations. Moreover, even when 
such topological phases are proposed to be realized in 
condensed matter materials, it is usually hard to directly 
image the wavefunction of, say, topologically protected 
bound states with current technology. Due to the out- 
standing controllability, artificial systems are promising 
alternative candidates for studying these novel phases, 
and there is a number of theoretical and experimental 
studies of topo logical phenomena using cold atoms and 
photons^HSl] Among them, quantum walks provide 
unique platform where any topological phase classified 
in one and two dimensions is realizable with simple mod- 
ifications of their protocols, as we will explain in this 
papeil^. In fact, one dimensional topological class pre- 
dicted to arise in Su-Schrieffer-Heeger model and Jackiw- 
Rebbi model has been already realized in the photonic 
architecture^. The key ingredient in the versatility of 
quantum walks is the controllability of the protocols. Be- 
cause protocols are something experimentalists choose to 
implement, it is possible to design the protocols in such 
a way to preserve or break certain kind of symmetries. 
From the following section, we study how such topologi- 
cal structure appears in discrete quantum walks. 
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IV. TOPOLOGICAL PHENOMENA IN 
QUANTUM WALKS 

A. Topological characterization of quantum walks 



In Section [II A| we have seen that the eigenvectors of 
quantum walks n(fe) illustrated in Fig. [2] have non-zero 
winding as k goes from ~ir to 7r. Such winding gives 
a topological characterization of quantum walks in the 
presence of sublattice (chiral) symmetry. 

Aas we noted in Section |II A the effective Hamilto- 
nian of the quantum walk possesses the sublattice (chiral) 
symmetry, which constrains the eigenvector of Hamilto- 
nian for each quasi-momentum n(fc) to lie on the plane 
perpendicular to the vector A. Under this sublattice 
(chiral) symmetry, n(fc) represents a map from the first 
Brillouin zone, which is a circle, to the equator of Bloch 
sphere, which is also a circle. Then quantum walk de- 
scribed by Eq. ^ possesses non-trivial winding of this 
map, where n(fc) goes around once the equator of Bloch 
sphere as k goes from — ir to ir. In the presence of 
the symmetry, the winding observed in Fig. [2] c) is ro- 
bust against small perturbations; one cannot change the 
winding number by small change of the Hamiltonian 
H c ff which preserves the sublattice (chiral) symmetry. 
One can intuitively check this robustness by trying to 
change the winding through the continuous deformations 
of n(fc). We call this winding number as topological num- 
ber due to their robustness against perturbations, and 
the topological number associated with iJ c ff is 1 whenever 
9 7^ 0, 27i^2l for our quantum walks. More generally, the 
winding of n(fc) around the equator can take any integer 
value, and different topological phases in this topologi- 
cal class are associated with different integers (winding 
numbers), and thus the topological classification is given 
by Z (a set of integers). 

It is important to note that this winding number is not 
topological in the absence of any symmetry constraint, in 
the sense that the winding number can be made zero by 
continuous change of n(/c). For example, one can shrink 
the loop of n(fc) into a point on a Bloch sphere. 

More generally, the concept of topological numbers 
(or topological invariants) are defined for a collection of 
Hamiltonian that represent band insulators with certain 
symmetries. A topological invariant is assigned to each 
band, and its value cannot change under the continuous 
deformations of Hamiltonian which preserve the symme- 
try There is one exception to this statement; when two 
bands mix with each other, the topological invariants can 
be changed in these two bands. We can flip this argu- 
ment and say that topological numbers can change their 
values only if two bands close their band gaps in the 
process. We argue in Section |IVC| that this property 
of topological numbers results in the creation of topo- 
logically protected bound states in the spatial boundary 
between regions that belong to two distinct topological 
phases. 

In the conventional quantum walk, the winding num- 



ber associated with the effective Hamiltonian is always 
I, and no other topological phase exist in this family of 
quantum walks. In the next section, Section |IVB| wc 
give yet another family of quantum walks that display 
two distinct topological phases with winding number 
and I in the presence of sublattice (chiral) symmetry. We 
illustrate how topological numbers can change as Hamil- 
tonian is changed in this example. 



B. Split step quantum walks 

In this section, we extend the conventional quantum 
walk by modifying the protocols and define so-called 
split-step quantum walks. This example illustrates how 
one can engineer topological phases in quantum walks 
through the active design of the protocols. 

Split-step quantum walks is a simple extension of con- 
ventional quantum walks which have one additional rota- 
tion and translation process (see Fig. [4] a)). The complete 
protocol is as follows; 

1. Rotation of the spin around y axis by angle 9i, 
corresponding to the operation R y (9i) = e~ l6lCr y/ 2 . 

2. Translation of spin up particle to the right, given 

by T t = Ej \j + 1)01 ® I tXt | + U)U !• Spin 

down particle stays in the original position. 

3. Second rotation of the spin around y axis by an- 
gle 92, corresponding to the operation R y {9 2 ) = 

e -i0 1 a y /2 

4. Translation of spin down particle to the left, given 

by n = Ej \j - 1)01 ® U)U I + 1 ® I t)<t |. Spin 
up particle stays in the original position. 

Thus, the evolution operator of one step is given by 
U(9 1 ,9 2 ) = T l R y (9 2 )T t R y {9 1 ). This split-step quantum 
walk is reduced to the conventional quantum walk de- 
fined in Section [II A| with 9 2 — 0. As before, we can find 
the effective Hamiltonian through U = e~ lHcti . The ef- 
fective Hamiltonian again takes the form Eq. ^ where 
the quasi-energy is 

cos E{k) = cos(6» 2 /2) cos(6»i /2) cos fc-sin(0i /2) sin(9 2 /2) , 

(11) 

and the eigenvector n(fc) is 



n x {k) 
n y (k) 
n z (k) 



cos(6» 2 /2)sin(6»i/2)sinfc 
sTn~B(fc) 

sin(0 2 /2) cos(0i/2) + cos(6» 2 /2) sin(0i/2) cos k 
sin E(k) 

- cos(0 2 /2) cos(6>i/2) sin k 
sin£;(fc) ~' 



It is straightforward to check that A(#i) = 
(cos(6>i/2), 0, sin(6*i/2)) is perpendicular to n(k) 
for all k. Therefore, the system possesses sublat- 
tice (chiral) symmetry with the symmetry operator 



a) Split-step quantum walk b) Topological phase diagram 




2.Translate ^ s P in U P 



3. Rotate ^ 02 




4.Translate^ Spin down 



Repeat 



FIG. 4: a) Protocol for a split-step quantum walk. Split-step quantum walk is defined for spin 1/2 particle in one dimensional 
lattice, and consists of four operations; 1. spin rotation around y axis by angle 61, given by R y (0i); 2. Translation of spin up to 
the right by one lattice, Tf ; 3. spin rotation around y axis by angle 62, given by R y (02); 4. Translation of spin down to the left 
by one lattice, T4.. The evolution operator of one step is given by U = TiR y (82)TfR y (9i). b) The topological phase diagram 
of split-step quantum walk for various rotation angles 9\ and 62- The phase is characterized by the winding number W, and 
split-step quantum walks realize either W = or W = 1. Since winding number is a topological number, it can change its 
value only when the band gap closes, and these gapless phases are denoted by solid black line (band gap closes at quasi-energy 
E = 0) and by red dotted line (band gap closes at quasi-energy E — tt). 



r(0i) = e -^ A (»i)V 2 . Notice that the this symmetry 
operation only depends on the first rotation angle Q\. 

The existence of sublattice (chiral) symmetry allows us 
to characterize the split-step quantum walk by the wind- 
ing number, denoted by W ', of n(fc) around the equator 
of Bloch sphere. Using the explicit expression for n(fc) in 
Eq. (p]), we find W = 1 if | tan(0 2 /2)/tan((9i/2)| < 1, 
and W = if | tan(0 2 /2)/tan(0i/2)| > 1. Thus the split- 
step quantum walk can realize different winding number 
for different rotation angles 81 and 82- We plot the phase 
diagram of split-step quantum walk in Fig.[4]b). 

For a given 6\, a set of Hamiltonian for varying val- 
ues of second rotation 82 has the same sublattice (chiral) 
symmetry, and thus {Hg 1 (S2)} defines quantum walks in 
the same topological class. Because the dependence of 
He 1 {&2) on 82 is continuous, topological nature of the 
winding number implies that the winding number is the 
same for wide range of the rotation angle 62 as phase 
diagram Fig. [4] b) shows. 

However, the winding number can change its value 
when the two bands close their gap. Such gapless points 
are given by the points | tan(6> 2 /2)/ tan(#i/2)| = 1, or 
62 = ±#1 , 2-7T ± 6*i denoted by solid and dotted lines in 



Fig. [1] b) . The mechanism behind the change of wind- 
ing numbers is the following. At the value of O2 where 
two bands close the gap, for example, 82 = Si, eigenvec- 



tor n(fc) in Eq. ( 12 ) becomes ill-defined at the quasi- 



momentum k corresponding to the degenerate points, 
because any superposition of spin up and down is the 
eigenstate of the Hamiltonian at that points. At this 
82 = Si, the winding number cannot be defined, and 
winding numbers of bands at H 6l (81 — e) and H Sl (81 + e) 
do not have to be the same. 

Note that this argument can be used to find the phase 
diagram in Fig. [4] b) without calculating the winding 
number at each point of the phase diagram. Since the 
winding number can change only when the bands close 
their gap, the phase transition between two topological 
phases is always gapless. Thus, in order to draw the 
phase diagram, it is only necessary to identify the gap 
closing points in the parameter space, and find the wind- 
ing number of the region bounded by the gapless phase 
lines. In quantum walks, the periodicity of quasi-energy 
allows the closing of the gap at either quasi-energy E = 
and E = tt as we saw in Fig.[2]b). These gap closing lines 
at quasi-energy E = and E = tt are denoted in Fig. [4] 
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E(k) E (k) E{k) 




K * t 




62+ 



FIG. 5: Illustration of the existence of bound states across the 
boundary of regions that belong to distinct topological phases. 
Here, we consider the inhomogeneous split-step quantum walk 
where the second rotation 82 (x) changes in space, and the 
winding number associated with the phase 82- = &2(x — > 
—00) is different from the winding number associated with 
the phase 82+ = &2(x — > 00). In both limit of x — > —00 and 
x — > 00, the bands are gapped. However, the winding number 
cannot change its value unless band gap closes, and thus it 
is expected that the band gap closes in the middle near the 
origin. States at E — that exist near x = must be localized 
since there is no state at this energy far into the left or into 
the right of the system. 



b) as solid and dotted lines, respectively. 

The topological structure of split-step quantum walk 
described above has a strong asymmetry between 8\ 
and 82, but this asymmetry is an artifact. One can 
shift the starting time of the quantum walk by uni- 
tary transformation and define an equivalent dynam- 
ics through U' = T^R y {ei)T i R y {e 2 ). It is straightfor- 
ward to show that this quantum walk has sublattice 
(chiral) symmetry given by T(0 2 ) — e ~" rA ( 02 )' <T / 2 with 
= (cos(#i/2), 0, sin(0i/2)). In this case, the quan- 
tum walks with constant 82 correspond to the Hamilto- 
nians in the same topological classes. 



C. Physical manifestations of 
topological band structure 

The non-trivial winding, or topological number, of the 
bands in quantum walks gives rise to a robust bound 
states at the boundary between two phases with different 
topological numbers. In the following, we first give an 
intuition behind the existence of such robust edge states. 

Here we consider the split-step quantum walks with 
inhomogeneous rotations in space, in order to create a 
boundary between quantum walks with different winding 
numbers. Here we take the first rotation 8\ to be homo- 
geneous in space and make the second rotation 8 2 space 
dependent, see Fig. [5j The effective Hamiltonian of this 
inhomogeneous quantum walk, having the homogeneous 
first rotation 6\, possesses the sublatice (chiral) symme- 
try given by Y{9 1 ) = e -«rA(0i)-o-/2_ while this statement 
is intuitively clear, it is instructive to explicitly show the 



existence of sublattice (chiral) symmetry, and we provide 
the proof in the Appendix [5] 

We take the second rotation angle to approach 62- for 
x — > — 00 and 82+ for x — > 00. We take the region in 
which the rotation angle changes from 2 - to # 2 + to be 
finite region around x — 0. Now we take 02—, 82+ to be 
such that winding number is for the split-step quantum 
walk with the rotation angle 8\, 82- and winding number 
is 1 for the walk with the rotation angle 8 1 ,d 2 +. Then 
the region near the origin represents the phase boundary 
between two distinct topological phases. 

In Fig. [5j we illustrate the local band structures in this 
inhomogeneous quantum walk. Strictly speaking, the 
system is spatially inhomogeneous, so quasi-momentum 
is no longer a good quantum number, and band struc- 
tures cannot be drawn. Yet, it is helpful to visualize band 
structures to understand what happens at the boundary. 
If the variation of rotation angle 82 is slow, then one can 
visualize the band structures at point xq to be the band 
structures of homogeneous quantum walk with rotation 
angles 8\ and 8 2 (xq) . Such description is certainly appli- 
cable in the limit x — > —00 and x — > 00. 

The band structures of both x —00 and x —¥ 00 rep- 
resent band insulators where two bands are separated by 
a band gap. By definition, the two bands are character- 
ized by winding number W — in x — > — 00 and W = 1 
in x — > 00. However, since the sublattice (chiral) sym- 
metry exists throughout the space, the winding numbers 
can only change by closing the gap across the boundary 
near x = 0. Therefore, the gap must close near x — 0, as 
illustrated in Fig. [5] 

This argument shows that there must be states in the 
gap (near E — 0) around the origin. Now because there 
is no state near E = in the limit x — > 00 and x — > 
—00 (this energy is in the gap of the bands), we can 
conclude these states near E = around the origin must 
be confined around the origin. Thus there is generically 
a bound state near E = at the boundary between two 
different topological phases. 

This prediction can be confirmed by running a sim- 
ple simulations of inhomogeneous quantum walks. The 
presence or absence of bound states can be confirmed by 
initializing the particle near the origin and running the 
quantum walk protocols. If there are bound states near 
the origin, there is generically an overlap between the ini- 
tial state and the bound state, and even after many steps 
of quantum walk, there remains a non-zero probability to 
find the particle near the origin. On the other hand, if 
there is no bound state, the particle quickly propagates 
away from the origin due to the ballistic propagation of 
quantum walks as described in Section [II C| 

In Fig. |6j we present the result of two inhomogeneous 
quantum walks. In Fig. [6] a), the boundary between 
two topologically distinct phases is created near the ori- 
gin with winding number W — as 1 — > —00 and 
W = 1 as x — > 00. Specifically, the uniform first ro- 
tation is 61 = —it/2 and second rotation 82- — 37r/4 and 
82+ = 7r/4. Here we considered a smooth variation of the 
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FIG. 6: Evolution of the spatially inhomogeneous split-step quantum walk. The initial spin of the particle is spin up, and 
particle is initialized at x = 0. a) The rotation angles of quantum walk are chosen such that the quantum walk corresponds to 
trivial topological phase with winding number W — as x — > — oo indicated as a white dot in the phase diagram and non-trivial 
topological phase with W = 1 as x — > oo indicated as a blue dot. Here we took the uniform first rotation 8\ = —n/2 and second 
rotation 82- = 37r/4 and 82+ = 7r/4 with $2 (x) = 5 (#2- + #2+) + 5(^2+ — O2-) tanh(x/3). After many steps of quantum walks, 
a large probability of finding a particle near the origin remains, indicating the presence of bound states, b) In this quantum 
walk, the phase of the two sides of the origin has the same winding number. The phase as x — > —00 is indicated by the white 
dot in the phase diagram and phase as x — > 00 is indicated by the red dot. Here we took 81 = —ir/2 and 62- = 37r/4 as before 
and 82+ = ll7r/8. In this case, the probability to find the walker near x — after many steps decays to 0, indicating the 
absence of a localized state at the boundary. 



second rotation given by 02(x) = 5 (#2- +^2+) + 5(^2+ — 
62-) tanh(a:/3), where the second rotation changes from 
62- to 6*2+ with the length scale of ~ 6 sites. The phases 
of quantum walks in the limit x — > —00 and x — > 00 are 
indicated on the phase diagram as the white and blue 
dot, respectively. In the simulation, the initial spin of 
the particle is spin up, and particle starts at x = 0. As 
we expect, a peak in the probability distribution appears 
even after 60 steps of the simulation, indicating the exis- 
tence of topological bound states. 

On the other hand, we studied the quantum walk in 
Fig. [6] b), where the system is characterized by W — 
throughout the space. Here we took 6\ = — n/2, 
62- = 3tt/A and 02+ = ll7r/8 again with the same func- 
tional dependence on x, 9 2 (x) — §(#2- + $2+) + |(#2+ — 
#2-) tanh(x/3). The phases of quantum walks in the 
limit x — y — 00 and x — ¥ 00 are indicated on the phase 
diagram as the white and red dot, respectively. With 
the same initial state of spin up, we implemented the 
quantum walk simulation, and the resulting probability 
distribution shows a fast decay of probability near the 
origin. After 60 steps, the probability near the origin de- 
cays close to zero, indicating the absence of bound states. 

For further details of the simulation, the interactive 



demonstration of inhomogeneous quantum walks is avail- 
able on the Mathematica demonstration website^^, 
where one can change the values of 9i and 62 and run 
the quantum walks. 

The topological class realized by the conventional 
quantum walk and by the split-step quantum walk is 
the same topological class as is proposed in Su-Schrieffer- 
Heeger model of polyacetylene and Jackiw-Rebbi model. 
This is the topological class in one dimension with sub- 
lattice (chiral) symmetry. In this respect, quantum walk 
acts as a quantum simulator of the topological phase. 
Since quantum walks are realizable with many different 
systems such as ions, photons and cold atoms, they allow 
the study of topological phases in a manner that is not 
possible in traditional condensed matter materials. The 
proposal to study topological phases and topologically 
protected bound states in split-step quantum walk was 
first proposed in [4], and later realized in experiments 
with photonic architecture^. The controllability of the 
experimental apparatus not only allow the direct imag- 
ing of the wavefunctions of topological bound states, but 
also allow to confirm the robustness of the bound states 
with parameter changes, the signature of topological ori- 
gin for the bound states. This topological bound state in 
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one dimensional system has not been directly observed 
in materials such as polyacetylene, and this photonic ar- 
chitecture provided the first experimental imaging of the 
bound states with topological origin in one dimension. 



D. Quantum walks with a reflecting boundary 

A special case of Hamiltonian with trivial topology is 
the vacuum, where the topological number associated 
with the system is zero. Thus, we can make a phase 
boundary by simply terminating the quantum walk with 
winding number W = 1. In this setup, the boundary ex- 
ists at the edge of the quantum walk, and it is expected 
that a bound state exists at this edge according to the 



general argument in Section IV C 



As we briefly mentioned in Section |II D[ such a setup 
is not unphysical, and intriguing realization of quantum 
walks with a reflecting edge has been suggested in [T5] 
through the understanding of a particle under the elec- 
tric field in discrete energy level structure as a quantum 
walk. The transition between different levels occurs as 
Landau-Zener process in this system, which corresponds 
to the translation operation in quantum walks. Thus, the 
ground state of the system acts as the reflecting bound- 
ary. They predicted the existence of a bound state near 
the ground state, which we can now understand as a 
topological bound state as a result of non-trivial topo- 
logical number of quantum walks. 

Here we consider the conventional quantum walk de- 
scribed by U = TR y (9) which extends from x = — oo up 
to x = 0. The quantum walk is terminated at x = 0. 
In order to conserve the particle number, we require the 
operation at the boundary to be unitary, i.e. the spin f 
particle needs to be reflected at the edge x = 0. Here we 
take the following operation at the edge x = 0; 

1. Rotation of the spin around y axis by angle 9, as 
in other sites, given by R y (9) = e - iea v/ 2 . 

2. Translation of the spin J, to site x = — 1. Spin f 
stays at x = and its spin is flipped to spin J, with 
phase accumulation e %v 

Explicitly, the operation at x = is 

U(x = 0) = T cdgc R y (6) 

Tcdgc = |-l)<0|®|i>U|+e^|0><0|®U)<t| 

In order to have a quantum walk in a topological class, 
it is crucial to have the sublattice (chiral) symmetry of 
the whole system. In particular, the sublattice (chiral) 
symmetry needs to be present for the evolution operator 
including the edge. If we denote the total evolution of 
the system with an edge as U x <o, then we require 



r-^or = u: 



,r<0 



(12) 



where r=e" i7rA CT/2 with A = [cos(0/2), 0, sin(0/2)]. This 
is a simple extension of the definition of sublattice (chi- 
ral) symmetry in Eq. (18]) to evolution operator. It is 



straightforward to check that the necessary and sufficient 
condition for the existence of sublattice (chiral) symme- 
try in U x <q is <p = 0, tt for the phase accumulated at the 
reflecting boundary. 



According to the general argument in Section IV C 
bound states exist near the boundary of x = 0. For this 
simple quantum walk, it is not very difficult to obtain 
the analytical solution of the bound state. The details 
of the derivation is given in the Appendix [C] We note 
that a straightforward extension of the derivation given in 
the Appendix should allow similar analytical solutions of 
bound states for the inhomogeneous split-step quantum 
walks. 

Here we describe the solution for the boundary condi- 
tion of tp = and rotation angle 9 — tt/2. The analytical 
solution shows that the bound state is at quasi-energy 
E = tt and the wavefunction takes the form 



l^f=°(-i)) 



l 



(-iy 



1 - V2 
1 



0<J 
(13) 



A = 



1 



log(V2-l) 

where Af is the normalization factor. Since the localiza- 
tion length A ss 1.1 is small, this bound state is tightly 
localized around x = 0. 

For the same boundary condition ip = but differ- 
ent rotation angle 9 = 5ir/2, the evolution operator 
U = TR y (9) is different from the one with the rota- 
tion angle 9 — 7r/2 by only a minus sign, i.e. TR y (9 — 
57r/2) = —TR y (9 = it/2). Thus, the same wavefunction 
in Eq. ( 13 ) is the bound state for this case as well, but 
now the quasi-energy of the bound state is E = due to 
the extra minus sign. 

On the other hand, for the reflecting boundary condi- 
tion with phase accumulation ip — tt with 9 = tt/2, the 
bound state exists at quasi-energy E = and wavefunc- 
tion is 



1 



1 - V2 
1 



o< j 



log(v/2 - 1) 

Generally, the quasi-energy of bound states for quan- 
tum walks with sublattice (chiral) symmetry is at E = 
or E = n. We will see in the next section, Section |TV E| 
that these energies represent special points where sublat- 
tice (chiral) symmetry provides the topological protec- 
tion of the states at these energies. 

For the special rotation angles of 9 = tt, it is possible to 
obtain the bound state solution by following the quantum 
walk for a few steps. Since this calculation is elementary, 
the existence of a bound state can be easily understood. 
This rotation angle corresponds to the rotation operation 



Ry(9 = 7T) 



Wy 

-1 

1 
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Thus, the rotation turns | f) ->• | |) and | i-) ->• — | t)- 
Let us take the phase accumulation upon reflection to be 
ip = 0. If wc initialize the particle at x = with spin 
down, the quantum walk follows the following evolution. 

|o>®|4.)4|o)®(-|t»3|o>®(-U>) 

Thus in this special case, |0) <£> | \) is an eigenstate of 
quantum walk operator localized at the edge. Since the 
state gains a minus sign after one step of the quantum 
walk, the quasi-energy of the bound state is tt, in accor- 
dance with the result obtained in the general analytical 
solution. 



E. Topological protection of the bound states: 
topological invariants 

The bound states resulting from topology studied in 
Section |IV C| and Section |IVD| are protected in a sense 
that they are robust against small changes in the quan- 
tum walk protocols. The logic is the following; since the 
bound states are the result of topological winding num- 
bers, and topological winding numbers cannot change 
their values unless the bands close their gaps, the bound 
states cannot disappear for a small change of rotation 
angles unless they are changed by a large amount such 
that the gap of the corresponding effective Hamiltonian 
closes. 

There is a more direct way to confirm such robustness 
by simply looking at the spectrum. In Section |IVD| we 
found that the energies of the topological bound states is 
always either E = or E = tt. Because the spectrum of 
the bulk (or the spectrum of the system without bound- 
aries) is gapped, the total spectrum of the system studied 
in the previous section look as in Fig. [7] a) and b) , where 
a single localized state sits at E = or E = tt and the 
bands of states away from E = or E = tt correspond 
to extended states in the bulk. 

Now we argue that the energy of the bound state sit- 
ting at the energy E — or E = tt cannot be changed by 
a continuous change of Hamiltonian which preserves the 
sublattice (chiral) symmetry. First of all, the presence of 
sublattice (chiral) symmetry implies that the states with 
energy E and —E have to come in pairs. In order to 
shift the energy of the bound state at E — by a small 
amount e, then, it is necessary to create two states at 
energy ±e. However, since a single state cannot be split 
into two, this is impossible, and the energy of the single 
state initially at E = 0(tt) is pinned at E = 0(ir). Ac- 
cording to this argument, the only way to remove such a 
zero (or tt) energy bound state is to change the Hamil- 
tonian until the bulk energy bands close the gap so that 
bulk states mix with the the boundary state. 

The structure of the spectrum illustrated in Fig. [7] a) 
and b) is generic. All the topological bound states exist 
at E — or E = tt. It is possible to assign topological 
numbers to these bound states, which give yet another 
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FIG. 7: General spectrum structure for inhomogeneous quan- 
tum walks with bound states, a), b) in Section IV D| we ana- 
lytically showed that there is a or tt energy bound state at 
the boundary of a conventional quantum walk. In a similar 
fashion, these bound states also appear in the inhomogeneous 
split-step quantum walks. These bound states energy are well 
separated from the extended, bulk states. Sublattice (chiral) 
symmetry requires that states with energy E and —E appear 
in pairs, and thus, a single and tt energy bound state cannot 
disappear unless bulk energy gap closes at or tt. c) Inho- 
mogeneous conventional quantum walks can possess and tt 
energy bound states, even though the winding numbers asso- 
ciated with the both sides of the boundary are zero. Such ex- 
istence of two flavors of topologically protected bound states 
is a unique feature of periodically driven systems and quan- 
tum walks. As is shown in Section |IVE| each bound state 
at or tt energy is associated with a topological number ±1. 
For a given system, a sum of these topological numbers Q° 
for E — and Q n for E = tt is a conserved quantity that 
cannot change its value unless the gap of the system closes. 



understanding of the topological protection of the bound 
states. These topological numbers are different from the 
winding numbers we assigned to the quantum walk band 
structures. These topological numbers are now assigned 
to the bound states themselves. The following considera- 
tion shows that the topological classification of the quan- 
tum walks with sublattice (chiral) symmetry is Z x Z, 
which means any integer numbers of E = and E = tt 
energy states are topologically protected. 

In the following, we consider general one dimensional 
systems with sublattice (chiral) symmetry, as in the case 
of quantum walks. Here we consider the bound states 
at energy 0. Analogous arguments can be applied to 
the bound states at tt. Suppose that there is Nq num- 
ber of degenerate bound states with energy 0. We label 
these states by with a' = 1---Nq. Let the sub- 

lattice (chiral) symmetry operator be given by F. Sub- 
lattice (chiral) symmetry implies that we have the anti- 
commutation relation between V and Hamiltonian, H, 
such that {T,H} = 0. As a consequence, T 2 commutes 
with H . When there is no conserved quantity associated 
with T 2 , it is possible to choose the phase of T such that 
r = 1. For example, in the case of quantum walks, we 
choose T ~ ie lA ' CTlr / 2 . Because of the sublattice (chiral) 
symmetry, we know that r|</?°,) is also an eigenstate of 
H with E = 0, so r represents a rotation within the 
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subspace of zero energy states, {| Then we can 
choose the basis of zero energy states such that they be- 
come eigenstates of T. We denote the zero energy states 
in this basis as {|"0a)} an d their eigenvalues under T as 
{Q° a }. Since r 2 = 1, Q° is either ±1. 

We now show that the sum of eigenvalues, Q° = 
SqQqi represents the topological invariant associated 
with zero energy bound states. We define the integer 
number Q° for zero energy bound states and Q v for tt 
energy bound states constructed in an analogous fash- 
ion, as 

Q° = E«i r i^> ( 14 ) 

oc 

Q T = E(C|r|C> (is) 

a 

where {IV'S)} are the tt energy bound states. 

In order to show that these quantities are indeed topo- 
logical invariants, we show that perturbations of the 
Hamiltonian which preserves the sublattice (chiral) sym- 
metry cannot mix two states both at the zero (tt) energy 
with the same eigenvalues of T. This implies that such 
perturbations do not lift the energies of these states away 
from or tt. Thus, Q° (Q^) number of bound states at 
energy E = 0(E = tt) cannot change under small defor- 
mations of the Hamiltonian. 

Let H' be a perturbation to the system such that 
{r, H'} = 0. Now we evaluate the matrix element of 
{r, H'} = in the 0(tt) energy states. The result is 

= (^ a \{T,H'}\^) 

(2W° a m%) for Q°=Q° 
~ I <V£TO£) - = for Q°^Q° 

This calculation shows that bound states with the same 
eigenvalues Q° a cannot mix. On the other hand, the same 
calculation does not give any constraint on the mixing of 
states with different eigenvalues Q Q a . Because one can 
break up any finite change of the Hamiltonian into suc- 
cessive changes of small perturbations, one can repeat 
this argument and show that the values Q° and Q v can- 
not change unless the bound states at and tt energies 
mix with the bulk states. 



F. Breaking of sublattice (chiral) symmetry 

Here we briefly comment on the perturbations of the 
Hamiltonian that breaks sublattice (chiral) symmetry. In 
the previous section, we give the argument that zero or 
tt energy bound states are protected as long as the bands 
in the spectrum do not close the gap. Now one can ask 
what happens if we consider the perturbations of Hamil- 
tonians that break sublattice (chiral) symmetry. Since 
no topological structure can be defined in the absence of 
sublattice (chiral) symmetry, there is no longer topolog- 
ical protection of the bound states. Yet, the statement 



that bound states cannot disappear until they mix with 
the bulk states remains true. Therefore, if we perturb 
the system that possesses topological bound states by 
adding small perturbations that break sublattice (chiral) 
symmetry, the existence of the bound states is still pro- 
tected by the existence of the gap. However, the bulk 
gap does not have to close to remove such a bound state; 
now the energy of the bound state can take any value in 
the absence of sublattice (chiral) symmetry and the state 
can be lifted away from or 7r energy. 



G. and n energy bound states ; topological 
phenomena unique to periodically driven systems 

Zero-energy bound state in one dimensional system 
with sublattice (chiral) symmetry has been known for al- 
most 30 years, and their existence was first predicted by 
Su-Schrieffer-Heeger model of polyacetylene and Jackiw- 
Rebbi modePSESl. On the other hand, we saw in previ- 
ous sections that quantum walks have two topologically 
protected bound states; and tt energy states. The ap- 
pearance of tt energy states is the result of the periodicity 
of quasi-energy. In return, quasi-energy is a property of 
periodically driven systems, and thus such appearance 
of two flavors of topologically protected bound states is 
a unique phenomenon to driven systems, which cannot 
occur in non-driven systems. 

In previous examples, only one of these states, and 
tt energy bound states, appears in a single system. From 
the argument given in Section |IV C| which of E = and 
E = tt bound states appears across the boundary of two 
topological phases is determined by whether the band 
gap closes at quasi-energy E = or E = tt. For ex- 
ample, consider the creation of a boundary between two 
different topological phases in the inhomogeneous split- 
step quantum walk, as we considered in Section [IV C| By 
choosing #2- and #2+ appropriately, one can either make 
the boundary between the two phases to be gapless at 
either E = or E = tt, as one can see from the phase 
diagram in Fig. [4] When the boundary closes the gap 
at E = 0(tt), the bound state appears at E = 0(tt) as is 
depicted in Fig. [7] a) and b) . 

Now consider the creation of the boundary between 
the phases with the same winding numbers by setting 
#i=0. This quantum walk whose evolution operator 
is given by U = X|i? a (# 2 )?t is nothing but the con- 
ventional quantum walk described in Section |II A| with 
initial time shifted. This time-shifted quantum walk re- 
alizes only a single phase with W = fp3. If we set O2- 
to be -2tt < 6 2 - < and 9 2 + to be < 8 2 + < 2tt, 
then the two phases corresponding to x — > —00 and x — > 
both have W — 0. From the point of view of the wind- 
ing topological number defined on the band structures, 
one expects no topologically protected bound states to 
be present. However, this inhomogeneous quantum walk 
possesses two topological bound states at quasi-energies 
E = and E = tt. 
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The existence of two bound states near the origin can 
be easily confirmed for the simple case of #2- = —n and 
6 2 + = 7r where the rotation 6* 2 _ is applied to all the sites 
x < and the rotation O2+ is applied to sites < x. 
These rotations act on the spins as 



t)->-U>,U>-Ht) 
t) -> U>,U> -> -I t) 



Now we consider the evolution of the particle after one 
step for a particle at site x = 1 with spin either up or 
down. This state evolves as 



|o>®|t)^|i)®|t)4|i)®U)^|o)®U) 
|o> ® 1 4.) 3 1 ) ® u> 4 |o) ® 1 1> 3 |o) ® 1 1) 

Thus it is clear that \ipE=o) = |0) <8> ^ (| t) + I D) is an 
eigenstate of the one-step evolution operator with quasi- 
energy and \i[>e=-k) = |0) <8> (| t) - ID) has quasi- 
energy tt. It is straightforward to check that any other 
states in this system has energy E = ±ir/2. Because this 
system possesses a single bound state at energy E = 
and E = tt, these states cannot be removed from these 
states until the gap of the bulk states closes, and thus 
these states must be present for any parameter values 
— 2ir < 82- < and < 82+ < 2tt. The general spectrum 
of such inhomogeneous quantum walks is illustrated in 
Fig. [7]c). As we noted in SectionllVEl these bound states 
are associated with topological numbers. In this walk, 
since 6± = 0, the sublattice (chiral) symmetry is given 
by the operator V = a x . Then the topological number 
associated with the bound state |V'_e=o) is nothing but 
the eigenvalue of T, so \ipE=o) has Q° = 1 and \tp E=n ) 
has Q 7 * = -1. 

Such coexistence of and tt energy bound state can 
also be observed in inhomogeneous split-step quantum 
walks, where the two phases on the left and on the right 
are separated by two gapless phases where one closes the 
gap at E = and the other at E = tt. 

The winding number is the comprehensive topological 
description of static Hamiltonians, but quantum walks 
are periodically driven systems. More completely, peri- 
odically driven systems should be described by the evolu- 
tion operator over one period, and the topological num- 
bers for such systems should be written in terms of the 
evolution operator and not in terms of the static effective 
Hamiltonian. Thus the topological classification of quan- 
tum walks is not given by Z as for the winding numbers, 
but in fact given by Z x Z as we have seen in the topo- 
logical invariants of bound states in Section IV E More 
detailed analysis of the difference of topological classi- 
fication between static systems and periodically driven 
systems is given in [36] . 



Two dimensional quantum walk 

1. Rotate 9l 2. Translate 
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FIG. 8: Protocol of two dimensional quantum walk consid- 
ered in Section [V A| The quantum walk is defined for a single 
spin 1/2 particle in two dimensional lattice. The protocol 
consists of 6 operations. In the figure, the spin-dependent 
translation is denoted by red arrows, where solid arrow is 
the movement for spin up and dotted arrow is for spin down. 
The protocols are 1. spin rotation around y axis by angle 
81; 2. spin-dependent translation where spin up is move to 
right and up by one lattice site, and spin down is moved to 
left and down; 3. spin rotation around y axis by angle 62; 
4. spin-dependent translation where spin up is moved to up, 
and spin down to down; 5. spin rotation around y axis by the 
same rotation angle as the first rotation 61; 6. spin-dependent 
translation where spin up is move to right and down to left. 
Each step of quantum walk takes a particle from even (odd) 
coordinate to even (odd) coordinate, so the lattice constant of 
the effective Hamiltonian is 2. Thus, the first Brillouin zone 
is -7r/2 < k x < tt/2 and —ir/2 < k y < tt/2. 



V. QUANTUM WALKS IN TWO DIMENSION 

A. effective Hamiltonian and Chern number 

In the previous sections, we illustrated the ideas of 
quantum walks and topological phases realized in these 
systems in the simplest setting, one dimensional quan- 
tum walks with two internal degrees of freedom. How- 
ever, the idea of topological phases is much more general, 
and it is possible to extend the quantum walk protocols 
to study different topological phases in different dimen- 
sions. In this section, we illustrate the idea by describing 
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the two dimensional quantum walks and demonstrating 
that this quantum walk realizes topological phases with 
Chern numbers, the phases that are responsible for inte- 
ger quantum Hall effects that we explained in Section [TTT] 
We consider the quantum walk of spin 1/2 particle on a 
square lattice. In the literature, quantum walks in dimen- 
sions larger than 1 are defined for larger number of in- 
ternal degrees of freedom, but the quantum walk defined 
here is simpler and easier to realize in experiments. The 
quantum walk consists of three rotations and three trans- 
lations, implemented in alternative fashion (see Fig. [8]); 

1. Rotation of the spin around y axis by angle 6\, 
given by R v (9i) = e'^y/ 2 . 

2. Translation of the spin j" one lattice to the right 
and up, and translation of the spin J, one lattice to 
the left and down. Explicitly, 



Ti 



y 

-i,y 



■ + l,y + l){x,y\® |t>(t 
l){x,y\®\l)(l\- 



3. Rotation of the spin around y axis by angle 02, 
given by R y {6 2 ) = e " i6 W2. 

4. Translation of the spin f by one lattice to up, and 
translation of the spin 4- by one lattice to down. 
Explicitly, 

72 = £ x ,>,y + i)<*,y|®|t)<tl 
+ \x, y -i){x,y\®\i){i\. 

5. Rotation of the spin around y axis by the same 
angle as the first rotation 9\, given by R y (9i) = 

g— i9l{Ty/2 

6. Translation of the spin t by one lattice to right, 
and translation of the spin 4- by one lattice to left. 
Explicitly, 

T» = E..> + i,y)(x,y\® I t)(t I 
+ \x-l,y)(x,y\»\i)(\.\. 

Note that in this quantum walk, the particle after one 
step of quantum walk moves from even (odd) coordi- 
nate sites to even (odd) coordinate sites as one can see 
from Fig. [HJ Thus the effect ive Hamiltonian of the quan- 
tum walk has the lattice constant equal to 2. There- 
fore, for translationally invariant quantum walks, the 
first Brillouin zone is given by — ir/2 < k x < n/2 and 
— 7r/2 < k y < ir/2. The evolution of the particle distri- 
bution in this walk can be studied in a similar fashion 
as the one dimensional analogue, and in particular, the 
asymptotic distribution is obtained through the formal- 



ism developed in Section II C 



In order to study the topological properties of this 
quantum walk, we consider the effective Hamiltonian of 
the quantum walk. As we detail in Appendix [D] the ef- 
fective Hamiltonian takes the form 



The spectrum £'(k)n(k) is determined by the equation 

cos(-E(k)) = cos(k x ) cos(k x + 2k y ) cos(#i) cos(# 2 /2) 

— sm(k x ) sin(k x + 2k y ) cos(#2/2) 

-ccs 2 (fc x )sin(0i)sin(0 2 /2) (17) 

The topological structure of two dimensional system 
appears in n(k) as in the case of one dimensional quan- 
tum walk. Since now we have two dimensional Brillouin 
zone, the function n(k) is a map from two dimensional 
torus to Bloch sphere, see Fig. [9] A small area on the 
torus is mapped to the small area on the Bloch sphere. 
If one maps the total area of the torus onto the Bloch 
sphere, the map necessarily wraps around the sphere an 
integer number of times due to the periodic boundary 
condition of the torus. Thus, if one calculates the total 
area covered by the map n(k), the value is 47m where 
n is an integer. This integer is so-called Chern num- 
ber, which is responsible for integer quantum Hall effect 
in two dimensional electron gas. Explicitly, the Chern 
number can be expressed in terms of n(k) as 



C = — / dk n • (d kT n x d k n) 
4tt./fbz 



(18) 



As opposed to the winding number of one dimensional 
quantum walk, this topological number does not rely on 
any symmetry of the system, and thus can exist in the 
absence of any symmetry. 

n • (<9 fc;c n x d ky n) 




First Brillouin zone 



Bloch sphere 



(16) 



FIG. 9: Illustration of Chern number of two dimensional 
systems with two bands. The eigenstate for each quasi- 
momentum k is a superposition of spin up and down, and 
can be represented as a point on Bloch sphere, given by n(k). 
Thus n(k) represents a map from the first Brillouin zone to 
Bloch sphere. In order to obtain a topological number for this 
system, we consider the area mapped by n(k) from the first 
Brillouin zone to Bloch sphere. Due to the periodic boundary 
condition of the first Brillouin zone, which is a torus, such 
map has to wrap around the Bloch sphere by integer number 
of times. This integer is what is called a Chern number, and 
represents the topological number associated with the system. 
The formal expression of the Chern number is then obtained 
by calculating the area covered by the gap n(k), which can 
be calculated in the way illustrated in this picture and results 
in the expression Eq. ( 18 1 
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J -7r/2 tt/2 

FIG. 10: a) Phase diagram of the two dimensional quantum walk. Each phase is characterized by a Chern number of a lower 
band with quasi-energy — -E(k), which can take values 0,±1 in this quantum walk. The Chern number can change only when 
the system crosses gapless phases, and the lines of gapless phases are indicated by the red and blue line in the diagram, b) 
physical manifestation of Chern numbers appears at the boundary between regions that belong to phases with different Chern 
numbers. Here we illustrate the inhomogeneous quantum walk, where the quantum walk in the central region (colored as red) 
corresponds to #i+ = 7n/6 and 82+ ~ 7ir/6 (red dot in a)), whereas in the other half (colored as white), the quantum walk 
corresponds to = 3n/2 and 82- = 3n/2 (white dot in a)). In the simulation, we took the periodic boundary condition 
for both x and y direction, and the system size is 100 x 100. c) Quasi-energy spectrum of the inhomogeneous quantum walks 
illustrated in b). The states colored as red are bulk states, and the states that go from the lower band to the upper band are 
the unidirectionally propagating edge states that are localized near the boundary of the two phases. States colored as blue are 
the states that run along the upper edge and those colored as green are the states that run along the lower edge, as illustrated 
in b). 



Conventionally, Chern number is associated with each 
band of Hamiltonian. The definition of Chern number 
above gives the Chern number of "lower" band with 
quasi-energy of —E(k) in Eq. ( 17 ), and the Chern number 
of upper band is given by simply — C so that the Chern 
numbers of all the bands sum to zero. More generally, if 
the wavefunction of a band at a given quasi-momentum 
k is given by |^(k)) = e tr ' u \ip(k)) , where |v(k)) is the 
periodic part of the Bloch wave function, then the Chern 
number associated with the band is given by the famous 
TKNN formulcP; 

C = J- f dk[d kx A ky -d ky A kx ] (19) 
L ~ Jfbz 

(A K ,A ky ) = (i(^(k)|^J^(k)),i^(k)|a fc J^(k))) 

This TKNN formula calculates the Berry phase of an 
electron as it goes around the first Brillouin zone. This 
expression of Chern number reduces to Eq. ( 18 1 in the 
case of systems with two bands. 

One can calculate the Chern numbers for the two di- 
mensional quantum walk described above for various val- 



ues of 0\ and 62- The phase diagram is plotted in Fig. 10 
a). A convenient way to obtain such phase diagram is 
to first obtain the lines of gapless phases. Since the 
topological number can only change its value across gap- 
less phase, it is only necessary to compute the Chern 
number at a single point of the gapped region and any 
gapped phase that is continuously connected with that 
point without crossing gapless phase must have the same 



topological number as that point. We detail the calcula- 
tion to identify gapless phases for this two dimensional 
quantum walk in the Appendix |E| 

A physical manifestation of topological phases with 
Chern number is, as is the case for one dimensional quan- 
tum walk, boundary states across the regions in which 
two different topological phases are realized. The exis- 
tence of such bound states can be understood according 
to the argument given in Section IV C the band struc- 



tures away from the boundary are gapped, but the band 
gap has to close near the boundary in order for the topo- 
logical number, Chern number, to change its value. Thus 
there is generically states in the gap of the bulk systems 
and these states are necessarily localized near the bound- 
ary. 

These bound states that appear in systems with non- 
zero Chern number are known to propagate in unidirec- 
tional fashion without any backscattering. It is possible 
to confirm the existence of such unidirectional edge states 
by considering inhomogeneous quantum walk where the 
particle at sites with < y evolves according to the two 
dimensional quantum walk with rotation angles 6\ + and 
82+, and the particle at sites with y < evolves accord- 
ing to the two dimensional quantum walk with rotation 
angles and #2-- Such inhomogeneous quantum walk 
is illustrated in Fig. 10 b). If the Chern number of the 
phases corresponding to 9\ + and 62+ and 6>i_ and #2- 
are different, unidirectional edge states are expected to 
appear along y = 0. 

Since these edge states exist, just like and 7r energy 
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states of one dimensional quantum walk, in the gap of the 
bulk states, it is easy to identify the existence of these 
states by numerically solving for the quasi-energy spec- 
trum. In Fig. 10 c), we provided the plot of quasi-energy 
spectrum for a torus geometry with periodic boundary 
condition on both x and y direction. The system size 
is taken to be 100 x 100. In the upper half of the sys- 
tem between < y < 50, we implemented the quantum 
walk with 8\ + — 7tt/6 and 82+ = 7tt/6, whereas in the 
lower half —50 < y < 0, the quantum walk corresponds 
to 6>!_ = 37r/2 and 82- = 3ir/2. Note that there are 
two boundary in this system, corresponding to the lower 
edge at y = and upper edge at y — 50. In the spec- 
trum, there are clearly two edge states, colored as red 
and green, which run along the upper and lower edge, 
respectively. These chiral edge modes are the signature 
of Chern numbers in two dimensional quantum walk. 



B. Unidirectionally propagating modes in quantum 
walks without Chern numbers 



In the case of one dimensional quantum walk, we found 
the existence of two bound states at quasi-energy and tt 
near the boundary of the phases with zero winding num- 
ber. This existence of two flavors of topologically pro- 
tected bound states represented a phenomenon unique 
to periodically driven systems that do not exist in static 
systems, and thus the existence is not captured by the 
winding number of the effective Hamiltonian. 

In a similar fashion, it is possible to have unidirection- 
ally propagating modes across the boundary of the re- 
gions where quantum walks in each region have no Chern 
numbers associated with the phases. In the following, we 
show that such chiral propagating modes exist for even 
simpler version of two dimensional quantum walk proto- 
cols. 

Here we consider the following simple two dimensional 
quantum walk with two rotations and two spin dependent 
translations, see Fig. [IT] 

1. Rotation of the spin around y axis by angle 9\, 
given by i? a (6q) = e -^i<y y /2_ 

2. Translation of the spin f one lattice to the right, 
and translation of the spin J. one lattice to the left. 
Explicitly, 

T i = J2x. y \x + ^,y)(x,y\ ® I t)(t I 
+ |&-l,y}<aM/|<8|4.}UI- 

3. Rotation of the spin around y axis by angle 82, 
given by R y {8 2 ) = e^ 2 " 7 */ 2 ; . 

4. Translation of the spin f one lattice to the up, and 
translation of the spin J. one lattice to the down. 
Explicitly, 

r a = Ex,vls.i/ + i>fofl®lt><tl 

+ \x,y-l)(x,y\<S\l)(\.\. 



Simple two dimensional quantum walk 

1. Rotate 01 2. Translate 



3. Rotate 02 



4. Translate 
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FIG. 11: Protocol for the simple two dimensional quantum 
walk studied in Section |V B| In this protocol, only four op- 
erations are applied during one step of quantum wallk. As 
before, the spin-dependent translation is indicated by red ar- 
rows, where solid arrow is for spin up and dotted arrow is for 
spin down. The explicit protocol is; 1. spin rotation around y 
axis by angle 0i; 2. spin-dependent translation where spin up 
is move to right, and spin down to left; 3. spin rotation around 
y axis by angle 82; 4. spin-dependent translation where spin 
up is moved to up, and spin down to down. 



The effective Hamiltonian of this quantum walk is 
again given by the form H e g = ^ k i£(k)n(k) -a® |k)(k|. 
The spectrum of this quantum walk is given 



cos(£'(k) = cos(k x + fc a )cos(6» 1 /2)cos(6» 2 /2) 



cos(fec — ky) sin(#i/2) sin(02/2) 



This quantum walk is described by Chern number zero 
phase everywhere, and the phase diagram is given in 
Fig. 12 a). All the gapless phases close their gap at both 
and 7T energy. Now consider the inhomogeneous quantum 
walks in this protocol, where the particle is controlled by 
the quantum walk protocol with rotation angles 6q + and 
82+ at sites < y, and the protocol at sites with y < 
is given by the rotation angles 81- and 2 -- If we choose 
the angles such that the two phases are separated by a 
single gapless phase, there are in fact two unidirectionally 
propagating modes at the boundary. 

This can be most easily confirmed for the spacial ro- 
tation angles 8\ + = and 82+ = tt and 8\- — ir and 
8<2- — by simply considering the evolution for spin up 
and down for a few steps. Near the boundary, the evolu- 
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FIG. 12: a) Phase diagram of simple two dimensional quantum walk studied in Section V B The Chern number of the quantum 
walk is everywhere zero. Yet, there are topologically protected unidirectionally propagating modes in the inhomogeneous 
quantum walk. The existence of such unidirectional edge states are shown through analytical calculations for special values 
of quantum walks in the text, b) quasi-energy spectrum of the inhomogeneous quantum walks, where the two quantum walk 
protocols corresponding to the two regions are indicated as orange and blue dot in a). The bulk states are colored as blue, and 
unidirectionally propagating states are colored as green and red. Green states propagate along the lower edge of the boundary 
and red states along the upper edge. For a given edge, say, lower edge, the edge states have non-zero energy winding as k x goes 
from — 7r to 7r, and thus, such edge states cannot be removed under the continuous change of quantum walk protocols unless 
the bulk gap closes, where the upper edge and lower edge are allowed to mix. 



tion is 

|j,o)®|t)4|j,o)®|t>4|j + i,o)®|t> 

4 |j + 1,0)® | |) 4 |. ? + 1,-1)® ||) 

\j, -l) ® U) 4 -i) ® |t) 4 -\j + 1, -i) ® 1 1) 
4-|j + i,-i>®|t)^-|j + i,o}®U) 

Thus we see that spin up states at site y — and spin 
down state at site y = — 1 both propagate to the right 
during the evolution. 

By Fourier transform in x coordinate, it is clear that 
the walk takes \k x ,y = Q)® \ f) -> e iK \k x , y = -1) ® | \) 
and \k x ,y = -1) ® | |) -> -e ik *\k x ,y = 0) ® | f). Thus 
we conclude there are two unidirectionally propagating 
modes 

|Vi) - -^-(\k x ,y = 0)®\t)+t\k x ,y = -l)<E>\l)) 
V i 

7T 

E\k%) = kx 2 

Itfe) = -j=(\k B ,y = 0)®\t)-i\k x ,y = -l)®\l)) 
E{k x ) = k x +* 

On the other hand, other states in the system evolve 



as, for I > 

>o) ® 1 1) 4 |j,z) ® 1 1) ^Ij + U) ® 1 1) 
4 |i + i,o®U) 4 \j + 1,1-1)® id 
4 |j + i,z-i> ®u> 4 |i,z-i)®|t) 
4-|j,z-i)®|t)4-|j,o®|t) 

Thus we conclude that the states 
5 ® 1 1) + i\j + I, I - 1) ® 1 1)) are eigenstates 
of the system with the flat quasi-energy E = —tt/2 and 
\ ® I t) - i\j + lj I - 1) ® 1 1)) are another set of 
eigenstates with quasi-energy E = tt/2. In a similar 
fashion, it is straightforward to show that the bulk states 
in / < — 1 have energies ±7r/2. 

Notice that the edge states obtained above has a non- 
zero winding in the energy direction as k x goes from — tt 
to 7T i.e. the states run from quasi-energy E = —n to 
E = 7r as k x goes from — 7r to tt. it is straightforward 
to convince oneself, by drawing the spectrum, that such 
energy winding cannot be removed under the continuous 
change of quantum walk protocols unless the bulk gap 
closes. Thus, the existence of these states is guaranteed 
for the inhomogeneous quantum walk which has 9\ + and 
82+ that are continuously connected with 9\ + = and 
62+ = tt and 61— and 62— that are continuously connected 
with = tt and # 2 - = 0, without crossing the gapless 
phase. 

As an example, we plot the quasi-energy spectrum 
of inhomogeneous quantum walk with # 1+ = tt/8 and 
O2+ = 77r/8 from < y < 50, and 61- = 7tt/8 and 
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02- = — 7r/4 from —50 < y < 0, and we again take peri- 
odic boundary condition for both x and y direction with 
system size 100 x 100. These two phases are connected 
with the limit 0\+ = and 62+ = tt, and #i_ = tt and 
62- = in a continuous fashion, as one can check. The 
spectrum of this system is plotted in Fig. 12 b), and 



one can see the existence of unidirectionally propagating 
modes on the two edges, colored as red and blue. One 
observes these edge modes in fact winds in the energy 
direction with non-trivial winding number. Such energy 
winding is in fact closely related to the phenomenon of 
Thouless pump®, and we refer the interested readers to 
the detailed analysis in |36| . 



VI. OTHER TOPOLOGICAL PHASES 

Different class of topological phases exist in various 
symmetries and dimensions, as is classified for non- 
interacting static 

Hamiltoniarpra As we have illus- 
trated the ideas through a few examples in this review ar- 
ticle, it is possible to realize any of the topological phases 
in 1 and 2 dimension through the variations of quantum 
walk protocols. We refer the interested readers to the 
article [4] for more complete analysis. 

Furthermore, we have given two examples of topolog- 
ical phenomena that are unique to periodically driven 
systems; and tt energy bound states in zero winding 
number phases and energy winding unidirectional edge 
states in zero Chern number phases^EH. These phenom- 
ena can also be extended to other classes of driven sys- 
tems. Recently it has been proposed that two flavors of 
Majorana Fermions can be realized at and tt energies 
in cold atoms^, in a similar fashion as and tt energy 
states of quantum walks. It is of great interests to study 
if other types of topological phenomena unique to driven 
systems can be realized in quantum walks. 



VII. CONCLUSION AND OPEN QUESTIONS 

In this review article, we studied topological phases 
appearing in quantum walks. After the introduction of 
quantum walks, we provided a thorough explanation of 
topological nature of quantum walks. We first associ- 
ated the quantum walks with winding numbers, and gave 
an intuitive argument for the existence of bound states 
across the boundary of the regions that belong to dif- 
ferent topological phases. We argued for the topologi- 
cal protection of bound states in two different point of 
view; one from the spectrum (gap) in the system and 
another from topological invariant associated with the 
bound states. These physics are illustrated through the 
explicit example of quantum walks in one and two di- 
mensions. We also explicitly demonstrated the existence 
of topological phenomena unique to periodically driven 
systems in one and two dimensions. 



There are many open questions that one can study in 
the field of quantum walks. For example, there is not 
yet an example of three dimensional quantum walks that 
realizes non-trivial topological phase. A simple exam- 
ple of such quantum walks are of interests, considering 
the exciteme nt in the field of three dimensional topologi- 
cal insulator d 21 l 22 l Moreover, three dimensional quantum 
walks with spin 1/2 has the possibility to realize Hopf- 
insulator first proposed by [17] . Since the realization of 
this topological phase is very difficult in condensed mat- 
ter materials, it is of great interests to explore the possi- 
bility of realizations in artificial systems such as quantum 
walks. 

Other open direction is provided by quantum walks 
in different geometries, such as hexagonal lattice. Since 
hexagonal lattice has three neighbors, study of hexagonal 
lattice quantum walk with spin 1 might provide interest- 
ing platform to explore unique quantum phenomena. 

Less concretely, the study of quantum walks with a few 
to many particles with strong correlation would be inter- 
esting to investigate. In particular, in the presence of 
frustrated hopping, there may be unique quantum phe- 
nomena such as the formation of spin-liquid phase. 

The author would like to thank a numerous insight- 
ful discussions with Mark S. Rudner, Erez Berg, Yutaka 
Shikano, and Eugene Demler. 



Appendix A: Asymptotic distribution of quantum 
walk 



In this section, we derive the intuitive result Eq. ( A2 ) 



which gives the asymptotic distribution of quantum 
walks. Here we consider a quantum walk initially pre- 
pared at site x = with initial state \s) such that 
|i) = \x = 0) ® I a). The evolution of a particle after 
each step is dictated by the effective Hamiltonian given 
by H cS = J dkE(k)n(k) -a® \k){k\ as in Eq. 

Since the particle is propagating under the non- 
interacting Hamiltonian H c g, it is natural to expect that 
the particle propagates in a ballistic fashion. Thus, the 
particle distribution has a well-defined form in terms of 
the variable X — x/N in the asymptotic limit. The dis- 
tribution of X, P(X), can be computed through 



P(X) = (S(x/N-X)) 

dse is(X-£/N) 



dse lsX / P -^/ N \ 



(Al) 



Thus we aim to obtain the expectation of so-called char- 
acteristic function e ~ lsx / N after A^th steps as N — > 00. 
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We first note that x = JV and thus 



-U&/N _ i , \V u\/,, , , ^,2(3 /N) 2 



1 + EH*)^Ij)01 + ( 



-isj/N _ 



2! 



Ij'Xj'I 



1 li>01 



E e_isj7JV i-?')0i 



dfc|ft + «/A0<*l 



Now we evaluate (e isx / N \ = 

^ e iH elt N e -isx/N e -iH ell N^^ Since ^ jg diagonal 
in quasi-momentum space, the evaluation is straightfor- 
ward. First of all, 

e iH eS N e -isx/N e -iH et sN 

dk \k + S/N)(k\ ® e ^E(k + s/N )n (k + s/Nya 
-iNE(k)n(k)-a 

dk\k + s/N){k\®e lsVkll[k) -' J 

In the last line, we took the expression in the low- 
est order in s/N. This can be confirmed through 
the expansion exp {NE(k + s/N)n(k + s/N) ■ a} = 
cos {NE(k + s/N)} + i sin {NE(k + s/N)} n(k + s/N) ■ 
a w cos (NE(k) + sv k ) + i sin (NE(k) + sv k ) n(k) ■ a to 
the lowest order in j^. 

It is now straightforward to evaluate the expecta- 
tion value of above expression in the initia l sta te \i) = 



J*^ -j=\k) <S)\s). Using the expression Eq. (|Al|), we ob- 



-TT V27T 

tain the final expression 



P(X) 



71 dkl 
2^2 



(I + (n{k) ■ a)) 6(v k - X) 



+ -(l-(n(k)-a))6(v k + X) 



Appendix B: Sublattice (chiral) symmetry of 
inhomogeneous quantum walk 

In this section, we give the explicit proof of sublattice 
(chiral) symmetry for inhomogeneous split-step quantum 
walks. In Eq. Q, we defined the sublattice (chiral) sym- 
metry in terms of Hamiltonian. This definition directly 
translates to the sublattice (chiral) symmetry on the evo- 
lution operator after one period U as 



r^c/r = [/t 



(Bl) 



We have shown in Section IV B that the (homogeneous) 
split-step quantum walk, U = T^R y (9 2 )T^R y (9i) pos- 
sesses the symmetry with the operator Tg 1 =e~' mA ' a / 2 
where A = (cos0i/2, 0, — sin#i/2). Here we write the 



subscript 9\ on the symmetry operator to emphasize the 
dependence on 9±. 

Here Tg 1 is a local operator and thus, we expect that 
the chiral symmetry is preserved even if 9 2 becomes in- 
homogeneous in space. 

In order to explicitly check this, we expand the evo- 
lution operator U = TiR y (9 2 )T^R y (9i) in the posi- 
tion basis, where we take the general case that 92 de- 
pends on space in an arbitrary fashion. For example, 
T T= EA 1 + °z)/2\x + 1) (ar| + (1 - <r z )/2\x) (x\. After the 
expansion, one obtains 



I' = ^ 1 -^Ry{9 2 {x + l))^^R v {9 1 )®\x + l){x\ 



+ l^R y (9 2 (x + l))!^^) ® \ X )( X + i| 
1 ^R y (9 2 (x + l)) 1 ^R y (9 1 ) 

+ l A^Ry(d2{x)) l -^R y {9 l )\ \x){x\ 



Now the sublattice (chiral) symmetry condition 
(r^ ) UTg 1 — W can be checked by comparing both 
sides of the equation for each position operators of the 
form \x)(x + a\ with a = —1, 0, 1. For example, compar- 
ing the both sides of the equation for the coefficients of 
|x + l)(x|, sublattice (chiral) symmetry requires that 



(TeJ-'^^RMx + l))- 



R~ 



2 



-Rl 



h(x + l)) 



2 

1 - (7 Z 



R y (9i)Tg 1 



It is straightforward to check that this equality indeed 
holds for any rotation 9\ and 9 2 {x + 1). Repeating such 
process for the coefficients of \x)(x + 1| and |x)(a;|, one 
confirms the existence of sublattice (chiral) symmetry for 
inhomogeneous quantum walks. 

Notice that the split step quantum walk with 9i = 
is effectively the conventional quantum walk described 
in Section |II A| The quantum walk becomes U = 
T^R y (9 2 )T^, which is unitarily related to the conven- 
tional quantum walk U con = TR y (9 2 ) — T^T^R y (9 2 ) by 
the shift of time. Therefore, the explanation above also 
provides the proof that the disordered conventional quan- 
tum walk where the rotation angle at each site is random 
possesses sublattice (chiral) symmetry. 



Appendix C: Analytic solution of the bound state for 
quantum walks with reflecting boundary condition 

In this section, we give the analytical solution of bound 
states for quantum walks with reflecting boundary con- 
dition studied in Section |IVD| Given the quantum walk 
with reflecting boundary condition whose evolution op- 



erator is U x <o in Section IV D we look for the bound 
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states near x = 0. Generally, such bound state can be 
written as 



IV*> = £(caU> + <*tlt»li) 

J<0 



(CI) 



The approach we take is to directly solve the eigenvalue 
problem 



U x < \i> b ) = e- iEb \i> b ) 



(C2) 



where E b is the quasi-energy of the bound state. Com- 
parison of the two sides of the equation above together 
with the normalizability of the bound state wavefunction 
allows the solution of the problem. 

The left-hand side of the equation gives 

U x <o\ip b ) 

= I>Jj-U> + E £i,tU' + l>t} + e^co,t|0,r) 
i<o j<-i 

= E ^+i^bU>+E^i.tb'.t>+e lv So,T|0,l> 



j<-i 



j<0 



where the tilde coefficients Cj^ t i are related to the origi- 
nal coefficients Cj^j, through the rotation R y {0) as 



c?4 



i? a (0) 



cos(0/2) -sin(0/2) \ 
sin((9/2) cos(6»/2) J 



c j4 



Comparison of the equation above with the right hand 
side of Eq. (CI) immediately gives 



Cj+14 3 ^ _1 



(C3) 



In matrix form, the first two equations can be rewritten 



as 



e-^ \ ( Cj - t \ _ / sin(0/2) cos(0/2) ^ +vr 

(0/2) -sin(0/2) ) { c u ) I e-^ II-.-. ) < 1 



C ^Ue ,£l cos^ + S55(W e 2l£b sin (^/ 2 ) |( '' + ■'- j , ur ■ , [ (C4i 
94 / V sin(0/2) cos(0/2) 




This last equation is a recursive equation that relates the 
coefficients at site j + 1 to site j. We denote the matrix 
that relates them as K, which is a matrix that appears 



on the right hand side of Eq. ( C4 1 



The behavior of wavefunction in the limit of x — > — oo 
is determined by the eigenvalues of the matrix K. They 
are given by 



K = 



v± = 



A± = 



A+ 
A_ 




>^/e~ 2iB b+e 2iE b-2cos(e) 
sin(e) 
1 



cos(.E 6 ) ± ^cos 2 (.E fc ) - cos 2 (6»/2) 
cos(0/2) 



where N± in the expression of v± are the normalization 
factors. 

Then, the amplitude of bound states wavefunction at 
site —j is given by 



c -j,t 
c-j'4 



Q 



X\_ 
X j _ 



Q- 1 



co,t 



i<-i 



where Q = ( v + v ) . One crucial observation is 
A + A_ = 1, and therefore |A + | < 1 when cos(Et) < 



and |A_| < 1 when cos(£'f ) ) > 0. The normalizabil- 
ity of the bound state wavefunction requires that the 
amplitude (co,^,co,;) T is proportional to v + (v_) when 
cos(_E;,) < 0(cos(i?h) > 0). No normalizable bound 
state wavefunction exists when |A + | = |A_| = 1, or 
cos 2 (£ b ) - cos 2 (0/2) < 0. 

Additional constraint on the amplitudes cq^,Cqx 



comes from the equation Eq. ( C3 1 , namely, e 



e I¥3 co,t- Solving these two conditions give us 



C04 



sm(0/2)e- iv 



',sm(E b ) T ^/cos 2 (E b ) -cos 2 (0/2) 



where minus sign is for u_ or when cos(£'b) > and plus 
sign is for v + or when cos(E b ) < 0. 

When ip = and < < 2tt, the solution exists at 
energy E b — it with (co,^,co,;) T oc v + . For ip = and 
2ir < 9 < 47r, the bound state energy is E b = and 
(co,t' c o,l) T oc v_. 

On the other hand, when ip = ir and < < 2tt, the 
bound state energy is E b = with (co^,co,;) T oc u_ , 
whereas when ip = tt and 2tt < 9 < 4-7T, the bound state 
energy is E b = n with (c ,f, c ,;) T cx v + . 

The bound state wave function found above de- 
cays on the length scale of ~ l/|log(A_)| = 

|log(l-| S in(9/2)|)-l°8(c°s(fl/2))| - TllUS tll<3 eXteilt ° f b ° UIld 
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state approaches oo as — >• 0, 2n. On the other hand, 
the bound state becomes most localized when 6 — tt. 



Appendix D: Spectrum of two dimensional quantum 

walk 



Here we give the details of how to compute the spec- 
trum of two dimensional quantum walk introduced in 
Section [V A| The method introduced here is general and 
can be easily extended to other protocols of quantum 
walks in, say, higher dimensions. 

The evolution operator of one step for the two di- 
mensional quantum walk can be written, in the quasi- 
momentum space, as 

U(k x ,k y ) = e^^e-^/V^e-^/ 2 

xe i(fcx + fey)crj e -i0i(T„/2 

Most general form of the effective Hamiltonian result- 
ing from spin 1/2 system is given by 

iJ cff (k) = E (k) + £(k)n(k) • a (Dl) 

This is true because a generator of two by two unitary 
matrix is 1 and Pauli matrices a. 

Then the spectrum can be identified by considering the 
trace of evolution operator because 

TV(Z7 (k x , k y )) = Tr (exp(-iff eff (k))) 
= Tr {e _iBa(k) (cos (E(k)) - i sin (E(k)) n(k) • a)\ 
= 2e-^ o(k) cos(S(k)) 

The explicit evaluation of the trace of U{k x , k y ) shows 
that E (k) = for all k and 

cos(.E'(k)) = cos(fc x ) cos(k x + 2k y ) cos(6*i) cos(02/2) 

— sin(fca;) sin(fc x + 2k y ) cos(8 2 /2) 

- cos 2 {k x ) sin(^i) sin(0 2 /2) 



Appendix E: gapless phase of two dimensional 
quantum walk 

In this section, we detail the calculation to obtain the 



line of gapless phase in the phase diagram of Fig. 10 a). 
The two bands of the system closes the gap when the two 
eigenvalues of Hamiltonian 



H oS (k) = £(k)n(k) • a 



(El) 



becomes degenerate. Since the eigenvalues of n(k) • a is 
±1, the quasi-energy of the states become degenerate if 
E(k) = —E(k), which happens if E(k) = or tt due to 
the periodicity of the quasi-energy. 

The gapless phase occurs at the values of Q\ and 02 such 
that the equation Eq. ( 17) has a solution of E(k) = or 



tt for some value of k x and k y . On the other hand, at 
E(k) — or 7r, the evolution operator U(k x ,k y ) takes 
the value 1 or — 1, respectively. 

A simple way to obtain such values of 9\ and 9 2 is to 
look at (1,1) component of the evolution operator axx = 
U{k x , k y )[l, 1]. First we study the lines of gapless phases 
where the gap is closed at the quasi-energy 0. Then the 
equation 1 = an gives the condition 



1 = e lK {1 



- cos(k x 



)cos(0 2 /2) 



sin(k x + 2k y ) ( 
2fc y )cos6»i cos(0 2 /2) 
cos k x sin B\ sin(02/2)} 



(E2) 



Note that 1 is the maximum magnitude that RHS of the 
above equation attains for any values of k x ,k y , 0i and 2 . 
Therefore, we can simply maximize the RHS in terms 
of the variables k x and k y or alternatively the variables 
k x + 2k y and k x . The argument separates few cases. 

If sin 9\ sin (02 /2) is non-zero, then cosfc^ = ±1. If 
we take the two orthogonal variables k\ = k x + 2k y and 
&2 = 2k x — k y , then the first two terms of the RHS of 
Eq. ( E2 ) is only a function of ki . By differentiating the 



absolute square of RHS in terms of k 2 , one can show that 
the extremum of this value is taken when cos k x = ±1. 

Suppose cos = 1. Then e iksc = 1 and the equation 
Eq. (E2) requires the first term sin(k x + 2k y ) cos($2/2) 
to be zero. If we take sin(fca; + 2k y ) to be zero, then 
cos(k x + 2k y ) = ±1, and the condition is reduced to 1 = 
003(6*1 + 6*2/2) for plus sign and —1 = cos(6'i — 6*2/2) 
for minus sign. Thus the gapless phase exist whenever 
6*i +6» 2 /2 = 2jm and 6*i - 6> 2 /2 = 2nn + tt. If 003(6-2/2) = 
0, then the equation is solved only at discrete points of 
Ox =tt/2 + 2im with 6> 2 = 'Sir + Arm or 6*i = 37r/2 + 2nn 
with 6*2 = 7r + Aim. These cases are already included in 
the condition above. 



Similarly consider cos k x 



T. Then e 



-1 



and the equation Eq. (E2) again requires the first term 



sin(k x +2k y ) cos(6* 2 /2) to be zero. If we take sm(k x +2k y ) 
to be zero, then cos^ + 2k y ) = ±1, and the condi- 
tion is reduced to —1 = 003(6*1 — 6*2/2) for plus sign 
and 1 = 003(6*1 + $2/2) for minus sign. Thus the gap- 
less phase exist whenever 6*i — 6* 2 /2 = tt + 2nn and 
^i + 6*2/2 = 2nn. Thus these conditions give exactly 
the same gapless phases as for cos k x = 1 . 

Now suppose that sin0i sin(02/2) is zero. A new condi- 
tion appears when sin(6* 2 /2) = 0. Then cos(6*2/2) = ±1- 
But now, we can satisfy the condition Eq. ( E2 ) by setting 
cos(fca; 



2ky 



0. Since this does not require for 6*i to 



take any particular value, we conclude that gap closes for 
the line of 6 2 = nn. 

A similar consideration for quasi-energy E = tt gives 
the condition that the gap closes at E = tt for 6*1+6*2/2 = 
27m + tt, 6*i — 02/2 = 27m and 02 = mr. These results 



lead to the gapless phases plotted in Fig. 10 a) 
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